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ABSTRACT 



Context. The yellow hypergiant stars (YHGs) are extremely luminous and massive objects whose general properties are poorly known. 
Only two of this kind of star show massive circumstellar envelopes, IRC +10420 and AFGL2343. 

Aims. We aim to study the chemistry of the circumstellar envelopes around these two sources, by comparison with well known AGB 
stars and protoplanetary nebulae. We also estimate the abundances of the observed molecular species. 

Methods. We have performed single-dish observations of different transitions for twelve molecular species. We have compared the 
ratio of the intensities of the molecular transitions and of the estimated abundances in AFGL2343 and IRC +10420 with those in 
O-rich and C-rich AGB stars and protoplanetary nebulae. 

Results. Both YHGs, AFGL 2343, and IRC + 10420, have been found to have an O-rich chemistry similar to that in O-rich AGB stars, 
though for AFGL 2343 the emission of most molecules compared with 13 CO lines is relatively weak. Clear differences with the other 
evolved sources appear when we compare the line intensity corrected for distance and the profile widths which are, respectively, very 
intense and very wide in YHGs. The abundances obtained for IRC +10420 agree with those found in AGB stars, but in general those 
found in AFGL 2343, except for l3 CO, are too low. This apparently low molecular abundance in AFGL 2343 could be due to the fact 
that these molecules are present only in an inner region of the shell where the mass is relatively low. 
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1. Introduction 

The Yellow hypergiants stars (YHGs) are among the most lu- 
minous (5.3 < logL[L ] < 5.9) and massive (Mj n ; t ~ 20 M ) 
stars (see, as general references, de Jager 1998, Jones et al. 1993, 
Humphreys 1991). These objects are thought to be post-red su- 
pergiants evolving bluewards in the HR diagram, but the details 
of such an evolution are still unknown. In at least a few of them, 
the stellar temperature is rapidly increasing. For instance, the 
spectral type of the hypergiant IRC +10420 has changed from 
F8Ia to A5Ia in just 20 yr (Klochkova et al. 1997, Oudmaijer 
et al. 1996, Oudmaijer 1998). Humphreys et al. (2002) showed, 
however, that the wind in this source is optically thick, suggest- 
ing that the apparent spectral type changes are due to variations 
in the wind rather than to interior evolution. 

Although it is thought that, during the red and yellow phases, 
these heavy stars eject as much as one half of their initial 
mass (e.g. Maeder & Meynet 1988, de Jager 1998), only two 
YHGs, IRC +10420 (= IRAS 19244+1 1 15) and AFGL 2343 (= 
IRAS 191 14+0002 = HD179821) are known to have very heavy 
circumstellar envelopes (CSEs). Those CSEs were detected in 
molecular line emission, dust-scattered light and IR emission 
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(see Hawkins et al. 1995, Meixner et al. 1999, Bujarrabal et al. 
1992, Humphreys et al. 1997, Neri et al. 1998, Bujarrabal et al. 
2001, Castro-Carrizo et al. 2001). 

Recent results from Castro-Carrizo et al. (2007) show that 
these circumstellar envelopes have several solar masses, very 
high expansion velocities (~ 35 kms _I ), and that such an en- 
vironment was formed in the last ~ 6000 yr. The envelope prop- 
erties there found are compatible with a mass loss driven by ra- 
diation pressure. However, in the case of IRC +10420 the ex- 
istence of infalling material suggests that other processes are 
also present in the inner parts of the envelope (Humphreys et 
al. 2002). Other properties of these molecular shells, like their 
chemical composition, have not been well studied. 

This paper is devoted to the study of the chemistry in the 
envelopes around IRC +10420 and AFGL 2343. We have ob- 
served rotational lines of several molecules and, in particular, we 
compare the line intensities and molecular abundances found in 
these objects with those of AGB CSEs and protoplanetary nebu- 
lae (PPNe). 

Both IRC +10420 and AFGL 2343 have very probably O- 
rich chemistry, in view of their OH maser emission (Likkel 1989, 
Reid et al. 1979; among evolved objects, only O-rich stars are 
found to emit in OH masers), silicate-rich circumstellar grains 
(e.g. Molster et al. 2002) and C-poor atmospheric composition 
(Klochkova et al. 1997, Thevenin et al. 2000). 
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Fig. 1. Observed spectra towards AFGL 2343. The intensity scale is given in units of main-beam brightness temperature (K). 




Castro-Carrizo et al. (2007) argued that AFGL 2343 was a 
YHG rather than a PPN as proposed by Josselin & Lebre (2001). 
As we will see, our molecular data confirm this, since the molec- 
ular line properties of these objects are quite similar and signif- 
icantly different from those usual in both PPNe and AGB stars. 
On the other hand, the nature of IRC +10420 is well known (see 
Jones et al. 1993, Oudmaijer et al. 1996). Our data will be used, 
in general, to determine the main chemical properties of the en- 
velopes surrounding IRC +10420 and AFGL 2343. 



2. Observations 

We have used the IRAM 30m telescope, at Pico Veleta (Spain), 
to observe mm-wave molecular lines in the yellow hypergiants 
IRC + 1 0420 and AFGL 2343 . The observed lines are: 

- In the 3mm band: C ls O 7=1-0, HCN 7=1-0, H 13 CN 7=1-0, 
SiO 7=2-1, 29 SiO 7=2-1, CS 7=2-1, CN JV=l-0, SiS 7=5- 
4, SO 7k=2 2 -1i, HC 3 N 7=10-9, HNC 7=1-0, and HCO + 
7=1-0. 

- In the 2mm band: SiO 7=3-2 and CS 7=3-2. 

- In the 1mm band: C ls O 7=2-1, HCN 7=3-2, SiO 7=5-4, 
29 SiO 7=5^1, CS 7=5^1, CN N=2-l, SiS 7=15-14, HNC 
7=3-2, and HCO + 7=3-2. 



In order to compare the line properties in these sources with 
those in other evolved circumstellar envelopes, we have also 
observed the O-rich AGB stars RX Boo and TX Cam, the C- 
rich AGB star IRC +10216, the C-rich PPN CRL2688, and the 
young C-rich PN NGC 7027. 

The coordinates of the star and some stellar properties of 
the program stars are summarized in Table 1. In IRC +10216 
and CRL2688, that are intense emitters, some other lines also 
appeared within the simultaneous receiver bands of the above 
lines. 

Our observations were performed in June 2000. SIS 3, 2, 
and 1mm receiver bands were used, often simultaneously. The 
receivers were tuned always in SSB mode. From frequent point- 
ing measurements, we expect errors ~ 3". The spatial resolution 
is 12-13" at 1.3 mm (taking into account the effects of pointing 
errors) and 22-26" at 3 mm. Beam efficiencies ranged between 
~ 0.75 at 3mm wavelength and ~ 0.5 at 1mm. 

The data presented here are calibrated in units of main beam 
Rayleigh-Jeans-equivalent antenna temperature, 7 m b, using the 
chopper-wheel method by observing hot (ambient) and cold 
loads (liquid nitrogen). The atmospheric conditions were good, 
with zenith opacities ; 0.2 at A=lmm. In addition, some results 
were compared with previous observations by Bachiller et al. 
(1997a,b), Bujarrabal et al. (1994a), and the catalog of standard 
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Fig. 2. Observed spectra towards IRC +10420. The intensity scale is given in units of main-beam brightness temperature (K). 

Table 1. Program stars. References: (H) Hipparcos, (1) Jones et al. (1993), (2) Bujarrabal et al. (1994a), (3) Skinner et al. (1998), 
(4) Jourdain de Mouizon et al. (1990). 
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line intensities for the 30m telescope (Mauersberger et al. 1989) 
in order to improve calibration. From the comparison between 
our different spectra and with respect with previous data, we ex- 
pect calibration uncertainties of about 20%. 

The detected spectra in our main sources, IRC +10420 and 
AFGL 2343, are shown in Figs. 1 and 2. Note that, in the case 
of CN N=2-l, only the two main groups of line components 
were observed, corresponding to 7=5/2-3/2 and 7=3/2-1/2. In 
Figs. A. 1-A.7 (in electronic version only) we present the spectra 
of the other observed sources: CRL 2688, NGC7027, RXBoo, 



TX Cam and IRC + 10216. As baseline profile, only straight lines 
were subtracted. 



3. Observational results 

We have obtained molecular line data of two YHGs, 
IRC +10420 and AFGL 2343, as well as for other sources, see 
Table 1. 

In Table 2 we show a summary of the new observational re- 
sults: peak intensity (K), rms noise and profile area (K km s _1 ). 
These parameters were in most cases calculated for a spectral 
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resolution of 1 MHz for the 3mm lines and of 2 MHz for the 
2mm and 1mm lines. The actual observations were performed 
with various spectrometers, whose resolutions were degraded if 
necessary to calculate these parameters. The upper limits are 5- 
cr level, calculated as follows: Integrated Area < 5 cr AV/ yN, 
being cr the rms noise for the observed spectrum, AV the equiv- 
alent velocity width, and N the number of channels within this 
velocity width. For the calculation of AV in the case of CN, we 
take the equivalent width of CN in IRC +10216 (Bachiller et al 
1997b) as the sum of the areas of all the hyperfine transitions 
divided by its T mb . This CN equivalent width of IRC +10216 is 
renormalized multiplying it by the ratio between AV for 13 CO 
7=2-1 in the star for which the limit is being calculated and AV 
for 13 COinIRC+10216. 

We have also considered the observations of various 
molecules from Bujarrabal et al. (1994a) and Bujarrabal, Alcolea 
& Planesas (1992), 13 CO 7=1-0 and 2-1 data from Bujarrabal 
et al. (2001), CN N=l-0 and 2-1 data from Bachiller et al. 
(1997b), and HCO + 7=1-0 data from Bachiller et al. (1997a), 
Sanchez Contreras el al. (1997), and Lucas & Guelin (1990). 
Note the interest of using data obtained with the same telescope, 
the IRAM 30m dish, in order to better interpret the line ratios. 
Taking into account data from the sources in these papers, the 
total sample of objects which data are used in our analysis, to- 
gether with the YHGs, is : 

O-rich AGB stars and red supergiant stars (RSGs): 
RXBoo, TXCam, RCnc, RSCnc, VYCMa, RCas, oCet, 
NML Cyg, W Hya, R Leo, VX Sgr, IK Tau, RT Vir, IRC -10529, 
IRC +1001 1, OH26.5+0.6, OH 44.8-2.3. 
C-rich AGB stars: LP And, UUAur, UCam, YCVn, SCep, 
VCyg, UXDra, VHya, CIT-6, CRL865, CRL3068, IRC- 
10236, IRC +10216, IRC +20370, IRC +30374, IRC +60144. 
O-rich PPNe: OH 17.7-2.0, OH 23 1.8+4.2, IRAS 17436+5003. 
C-rich PPNe: CRL618, CRL2688, NGC7027, 
IRAS 07134+1005, IRAS 19500-1709. 

Note that the RSGs VYCMa, NML Cyg, and VXSgr are 
also included, but very few data on them are available. From 
the point of view of the line intensity ratio comparison, and also 
due to the lack of data for these objects, both RSGs and O-rich 
AGBs behave in a similar way. Also, NGC 7027, a C-rich young 
PN, is considered in the following discussion as a PPN. 

3. 1 . Line intensity comparison 

In order to compare the molecular emission in YHGs with that 
of other CSEs, we have represented in Figs. 3 to 5 integrated 
intensities of several pairs of relevant lines. Open symbols rep- 
resent O-rich objects, and filled symbols C-rich objects. Squares 
and triangles respectively represent detections and limits for the 
AGB stars (AGBs). The RSGs are plotted as O-rich AGBs with 
thicker symbols. Four-peak asterisks and three-peak asterisks 
represent, respectively, detections and limits for PPNe. These 
polygons with three vertices have one of them pointing leftwards 
or downwards, indicating for which transition we have a limit. 
Objects with limits for the two represented lines are not plotted. 
For the YHGs, a five-peak open asterisk represents AFGL2343, 
a six-peak asterisk is used for IRC +10420, and the limits are 
represented by an arrow from the center of the symbol. To make 
the comparison easier with the other objects, straight lines (here- 
after average lines) represent the points with a ratio equal to the 
geometric average of the profile area ratios for the two YHGs. 

In Fig. 3 we show ratios between integrated line intensities 
for 13 CO and other abundant molecules. In Fig. 4 and 5 we show 



line ratios between pairs of molecules other than 13 CO. In Fig. 3 
we can see that, in general, AFGL2343 and IRC +10420 lay in 
significantly different regions of the diagrams. For IRC +10420 
these line ratios are comparable to these typical in O-rich AGB 
stars. For AFGL 2343, however, they are weaker by factors rang- 
ing between ~ 2 and 10. On the other hand, in Figs. 4 and 5 both 
YHGs fall very close to the average line and to the region of 
the diagrams occupied by O-rich AGB stars, and their line ratios 
are clearly different from those of C-rich stars and PPNe. The 
situation of AFGL 2343 in Fig. 3 could be due to an overabun- 
dance in 13 CO, an underabundance of the rest of the molecules, 
or to a small mass of the emitting regions for these molecules. 
Note that, for optically thin lines, the ratio of the profile areas 
depends on the abundance of the molecule and the mass of the 
emitting region. We can conclude from these diagrams that the 
YHGs show, as expected, O-rich chemistry similar to that of the 
O-rich AGB stars. 

The relative weakness of most molecular lines with respect 
to 13 CO in AFGL 2343 is particularly noticeable for 28 SiO, 
29 SiO, HCN, CN and HNC. The 28 SiO/ 13 CO intensity ratio in 
AFGL 2343 is similar to that of PPNe, and lower than those ra- 
tios obtained for AGBs, O-rich or C-rich. We also confirm the 
previously-reported low SiO intensity in PPNe. 

The HCN intensity of AFGL 2343, relative to 13 CO, lays un- 
der all O-rich AGBs. The HNC/ 13 CO and CN/ 13 CO line inten- 
sity ratios in this source are the lowest among all sources in our 
sample in which HCN or CN is detected. 

To go a step further in the comparison between the molecu- 
lar line intensity in these kinds of stars, we will now use profile 
areas corrected for distance. In general, for the same intrinsic lu- 
minosity, the observed brightness intensity of sources at different 
distances from the observer follows a ~ D 2 law, where D is the 
distance. Multiplying the intensities obtained by D 2 , in kpc, we 
set, in some way, all the stars at the same distance of 1 kpc. The 
distances for the YHGs were taken from de Jager (1998) in the 
case of IRC +10420, and from Hipparcos parallax measurements 
for AFGL 2343, which, despite of its low accuracy, is compatible 
with the luminosity of a YHG. For the rest of stars the distances 
are give in Table 1 or were taken from Bujarrabal et al. (1994a). 

The result of the distance correction is presented in Fig. 6, 
for the following representative transitions for each chemistry 
type: SiO 7=2-1 (O-rich) and 7=3-2, HCN 7=1-0 and HNC 
7=1-0 (C-rich), all versus 13 CO. We found, in order of increas- 
ing line intensities, first the AGBs, then PPNe, and finally the 
YHGs. This difference is higher in the case of SiO, an O-bearing 
molecule, than for HCN and HNC, which are C-bearing. The 
reason for this behavior is that, as we have seen, the YHGs show 
O-rich chemistry. 

We note that, although the use of distance corrected areas 
shows a difference between YHGs and the other objects, it de- 
pends directly on the quality of the distance determination. If, as 
claimed by Josselin & Lebre (2001), the distance for AFGL 2343 
is smaller than that used here, it would lay in the AGB or PPN 
region of the diagrams in Fig. 6. 

3.2. Line widths 

The equivalent velocity width, AV, provides an independent 
comparison between YHGs, AGBs and PPNe. We have used as 
reference line 13 CO 7=1-0, for which we measured AV for all 
objects. 

The mean values found are: 
AGBs: AV = (22±9)kms- 1 
PPNe: AV = (29 ± 9) km s" 1 
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Table 2. Peak T m b, rms noise, and integrated area from our observations. The resolution is 2 MHz for the 1mm and 2mm data, and 
1MHz for the 3mm data. *: 2MHz resolution. **: 4MHz resolution. ~: Tentative detection. 13 CO data from Bujarrabal et al. (2001). 
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In Fig. 7 we show a histogram of the obtained equivalent 
line width. This statistical approach to AV in these evolved ob- 
jects signalizes a major difference between YHGs and AGBs and 
PPNe. Note that this comparison is distance independent. This 
result gives us another reason to claim that AFGL 2343 is prob- 
ably a YHG. 

Due to the lack of data for the RSGs we can not obtain 
an estimate of the equivalent width using this method. From 
12 CO 7=2-1 line data (Cernicharo et al. 1997) a main value of 
55 km s~' is found for RSGs, similar to that of YHGs. This is a 
confirmation of the evolutionary connection between these ob- 
jects (see e.g. de Jager 1998). 

Note that, for many PPNe, the lines show a core (region with 
low velocity) and wings (high velocity). The emission in the pro- 
file wings comes from regions with a high expansion velocity, 
even larger than 100 kms -1 (see e.g. Bujarrabal et al. 2001), 
but weaker than that of the profile core. The equivalent width is 



therefore dominated by the core region and it is not very large. In 
any case, the shape of the line profiles in PPNe is very different 
from that of AFGL 2343 and IRC+ 10420, which do not show 
wings (see Figs. 1 & 2). 

4. Abundances 

4.7. Abundance estimate method 

The formulation used to calculate the abundances from the mea- 
sured intensities is similar to that used by Bujarrabal et al. 
(2001). We suppose that the populations of all the rotational lev- 
els can be described by a single temperature T rot . We do not as- 
sume optically thin emission. The abundance of a molecule can 
be estimated from the equation: 



X = In 



S(T M ) Os 



8ttv 3 Q e E,/kT " 
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Fig. 3. Profile area diagrams showing some representative line transitions versus 13 CO 7=2-1. 



M 



(1) 



where X is the fractional abundance, Q is the partition function, 
Qb is the HPBW beam size, Q.$ is the size of the source, AV 
is the equivalent line width of the emission (in cm s _1 ), g u is the 
degeneration of the upper state, A is the Einstein coefficient, £7 is 
the energy of the lower state of the transition, and M is the mass 
of the emitting gas. For linear molecules, E\ = hB„Ji(Ji + 1), 
where B a is the rotational constant. Note that, in general, the 
source is much smaller than the beam and the lines are optically 
thin, therefore the abundance does not depend on Q5 . 

In Eq. 1 the logarithm term appears to account for optically 
thick emission. However this approximation is only useful for 
moderate opacities, r. For values of the optical depth larger than 



~ 2 , the results become too strongly affected by observational 
and model uncertainties. In other words, for high r, the emission 
will only come from the external layers of the envelope, since we 
are unable to detect the emission from the inner zones, and the 
derived abundance will just be a lower limit. So, when our cal- 
culations indicate an optical depth larger than 2, we just give the 
value corresponding to r = 2 as the lower limit to the abundance 
(this happens for all the transitions of SiO, 29 SiO, and HCN in 
IRC +10420) . We also note that this expression becomes merely 
approximate for t>1 , since we use mean values of r within the 
profile to avoid the integral in the observed velocity. 

Due to the low values of r rot we found (see below), we de- 
cided to numerically calculate the partition function with high 
accuracy, rather than to use approximations. The sum was pur- 
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Fig. 4. Profile area diagrams for some pairs of transitions of O-rich molecules versus some of C-rich ones. 



sued until the next term had reached a very low tolerance level. 
In the case of SO, Q was calculated using the energy levels given 
by Omontetal. (1993). 

CN has fine and hyperfine structure (see Figs. 1 & 2). In or- 
der to simplify the calculation of Q, we reduce this case to a two 
AMevel system by recombination of their Einstein coefficients. 



A u l - — y^ l gu^ua',la- (2) 

^ a8u a', a 

This procedure can be safely used only for optically thin 
emission, where the relative intensity of each line is proportional 
to the Einstein coefficient, A. 

The radius of the molecule-rich nebulae around AFGL 2343 
and IRC +10420 are taken from results from interferometric CO 
maps by Castro-Carrizo et al. (2007). We will assume that the 
13 CO emission extends as far as 12 CO, while the emission from 
the rest of the molecules is supposed to come from the densest 
part of the shell. In fact, in the case of IRC +10420, the SiO 
emission is known to appear in this dense region (Castro-Carrizo 
et al. 2001). The mass of the emitting region is also derived from 
the results obtained by Castro-Carrizo et al. (2007). These data 
are summarized in table 3. 

We estimate r rot from the line intensity ratio for the 
molecules in which we have observed more than one transition. 
Sometimes, the value of T rot found in this way was slightly lower 
than the observed brightness temperature, after correcting the 
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Fig. 6. Distance-corrected profile areas diagram for some repre- 
sentative lines versus 13 CO 7=2-1 
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Table 3. Assumed radii and mass of the emitting regions, based 
on the 12 CO data from Castro-Carrizo et al. (2007) 



Molecule 


«m(cm) 


ffout(cm) 


M(M ) 


ffi,(cm) 


S,„,(cm) 


M(M B ) 




AFGL2343 


IRC + 10420 


"CO 


no 15 


5 10" 


4.02 


310"" 


5.210" 


0.7 


Others 


1 10 15 


2.5 10 17 


3.31 


3 10" 


1.24 10 17 


0.24 



o Ll 



20 _ 1 40 

Av (km s ) 



Fig. 7. Histogram of the profile equivalent line width. Upper 
panel shows the comparison between AGBs and YHGs. Lower 
panel refers to PPNe versus YHGs 



observed T mb for the dilution factor, Q. s /Cl B . In these cases, we 
imposed T mt to be equal to the brightness temperature. 

For the molecules with only one observed transition or with 
only one detection, we took an average rotational temperature, 
which was found to be ~ 10 K for AFGL2343 and ~ 1 1 K for 
IRC +10420. In some cases, where we have a limit for the tran- 
sition with the highest 7, we obtained that the upper limit to the 
abundance from the non-detection is lower than the abundance 
derived from the detected line, using the mean T TO i. In this case, 
we decrease the rotational temperature until both values of the 
abundances derived are equal, being an upper limit for T mt . 

4.2. Abundance results 

We present the results of our abundance calculations in Table 4. 
We also show the derived rotational temperatures. The tempera- 
tures are systematically higher for the molecules that show high 
optical depths, as expected when the excitation temperatures are 
lower than the kinetic temperature. 

We deduce the mean densities for the CO emitting region 
from the data in Table 3: n ~ 5 x 10 3 cm 3 for AFGL2343, and 
n ~ 8 x 1 2 cirT 3 for IRC + 1 0420. The mean densities derived for 
the emitting regions of the other molecules are n ~ 4 x 10 4 crrT 3 
for AFGL2343 and n ~ 2 x 10 4 cm 3 for IRC +10420. These 
relatively low densities are compatible with the idea that most 
lines, except in particular those of CO, should be underexcited. 

AFGL2343 shows a general underabundance in all the 
molecules but 13 CO. The abundances of molecules other than 
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Table 4. Mean abundances and rotational temperatures derived 
for AFGL2343 and IRC +10420. forced to be equal to the 
average r rot value. 2 : T mt imposed to be equal to brightness tem- 
perature. 3 : from intensity upper limit; see text for details. 



Line 


<X)(AFGL2343) 


r rot (K) 


<X)(IRC + 10420) 


r I0t (K) 


13 CO 


2.6 10- 5 


9 


4.7 IO" 5 


10 


C 18 


1.5 icr 6 


20 


<6.6 10~ 6 


n 1 


CN 


3.6 10- 9 


6 


3.8 IO" 7 


6 


cs 


8.3 10- 9 


10 


1.9 10" 7 


9 


H 13 CN 






2.9 IO" 7 


ll 1 


HC 3 N 


5.5 10- 9 


10 1 


5.7 10- 8 


ll 1 


HCN 


1.1 io- 8 


7 


>2.2 IO" 6 


36 2 


HNC 


2.0 IO" 9 


<8 3 


1.6 IO" 7 


18 


HCO + 


2.0 IO" 9 


<7 3 


2.7 IO" 8 


<6 3 


SiO 


5.4 IO" 8 


11 


>1.310" 5 


79 2 


SiS 


8.5 IO" 9 


10 1 


1.5 IO" 7 


ll 1 


so 






1.6 IO 6 


11 


29 SiO 


6.9 IO" 9 


8 


>i.oio- 6 


15 2 



CO are in general more than ten times higher in IRC +10420 
than in AFGL2343. AFGL2343 shows an abundance of 13 CO 
comparable to that typical of AGB stars, while this abundance in 
IRC +10420 is higher by a factor 2. 

The emission of molecules in AFGL2343, apart from 13 CO, 
could come from a region smaller than that assumed here 
(Sect. 4.1), containing significantly less mass. This would re- 
sult in an increase in the derived molecular abundances (in that 
region) for AFGL2343. In fact the emission of such species 
in AGB envelopes usually comes from regions at a distance 
from the star of several 10 16 cm, although the SiO emission 
from IRC +10420 is confirmed to appear at a distance of about 
10 17 cm. We note the difference in the assumed radii of the dense 
region for both stars (see Table 3), being that of AFGL2343 
twice that of IRC +10420. 

4.3. Comparison with previously published results. 

With the results in Table 4, we are able to compare the 
abundances here calculated with other published results for 
YHGs, like those obtained in Bujarrabal et al. (1994a) and 
Bachiller et al. (1997b). In these papers the distance assumed for 
IRC +10420 is 3.4 kpc and for AFGL 2343 is 6 kpc. To compare 
the results of the abundance calculations first we must correct 
for the effects of the different assumed distances. 

The average difference is ~20 for IRC +10420 and around 
~300 for AFGL 2343. In order to find out the origin of these dis- 
crepancies we must focus on the assumptions made in the differ- 
ent calculations. Both Bujarrabal et al. (1994a) and Bachiller et 
al. (1997b) used the same method, assuming optically thin emis- 
sion. In fact, Eq. 1, here used to calculate abundances, would be- 
come that used by Bujarrabal et al. (1994a,b) and Bachiller et al. 
(1997b) if we assume, as in those papers, optically thin emission 
and a negligible background temperature (7b g ~ 0). For those 
cases with optically thin emission, the differences in the calcula- 
tions can only be due to the assumed values of the emitting mass 
and of r rot . For optically thick cases, the differences also depend 
on the source size and opacity. 

In Table 5 we can see the radii used and the masses ob- 
tained for the emitting regions by Bujarrabal et al. (1994a) and 
Bachiller et al. (1997b). The radii assumed by these authors for 
regions emitting in molecules other than CO in YHGs are signif- 
icantly lower than the values obtained from PdB maps (Sec. 3), 



leading to smaller masses and, therefore, to higher abundances. 
Indeed, the discrepancies found in the abundances are almost 
fully explained by these differences in the masses adopted for 
the emitting region. 



Table 5. Radii and mass of the emitting gas for each molecule 
used in Bujarrabal et al. (1994a) and in Bachiller et al. (1997b). 
*: Shell, R m = 3.1 10 16 cm. 



Molecule 


R out (cm) 


M(M ) 


R out (cm) 


M(M ) 




AFGL 2343 


IRC+ 10420 


IJ CO 


10" 


0.2 


10" 


0.06 


SiO 


210 15 


4.2 10~ 3 


210 15 


1.2 10~ 3 


CN 






1.1 10 17 * 


0.05 


Other 


10 16 


0.02 


10 16 


610- 3 



Another important factor in the abundance estimate is the 
rotational temperature. In Bujarrabal et al. (1994a,b), a T I0t of 
20 K is assumed. The temperatures we found here are in most 
of the cases lower than 20 K. Bachiller et al. (1997b), however, 
estimate the rotational temperature from the ratio between CN 
N-2-l and N- 1-0 emission, and it is similar to the one obtained 
here for this molecular transition. 

In some cases, the ratio between the abundances calculated 
here and those from the cited authors is lower than the ratio of 
the masses. This is explained by the difference in T rot . Lower 
rotational temperatures lead to higher abundances in the case of 
transitions with a high J, like in the case of CS J =3-2 and J = 
5-4, which were the transitions used to calculate the abundance 
in Bujarrabal et al. (1994a). 

4.4. Comparison with abundances in AGB stars. 

The relative abundances, X, determined here (see Table 4) must 
be compared with the standard abundances found in circumstel- 
lar envelopes (e.g. Bujarrabal et al. 1994a), in order to under- 
stand the differences between our YHGs and those. 

The values of X found for IRC +10420 are in full agreement 
with those supposed for the O-rich CSEs, except for the high 
abundances obtained for 13 CO and HCN. In the case of 13 CO, 
this overabundance is a factor 2. Note that the value of Xhcn is 
affected by the opacity, as said in Sect. 4.1, therefore low abun- 
dances were found by Bujarrabal et al. (1994), who assumed low 
optical depths. 

In AFGL 2343, except for 13 CO, the abundances are much 
lower (around 70 times on average) than the standard values for 
the CSEs around AGBs. As we mentioned in Sect. 4.2., such low 
abundances found in AFGL 2343 could be related to the extent 
of the circumstellar layer emitting in these molecular lines. In 
AFGL 2343, the region we assumed to be molecule rich is as 
extended as 2.5 10 17 cm, but there is no observational confir- 
mation of such a large extent. If the region at which molecules 
other than CO are abundant in this source is comparable to that 
of IRC +10420, then the emitting mass would be much smaller 
than the mass assumed in our abundance calculation (~0.05M o , 
see mass estimates by Castro-Carrizo et al. 2007). This would 
yield molecular abundances in this region comparable to those 
usually found in molecule rich CSEs. 

5. Conclusions 

We observed several molecular transitions of a wide variety 
of species (see Sect. 2) in the yellow hypergiants (YHGs) 
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AFGL2343 and IRC +10420 as well as in several AGB stars 
and PPNe. Previously published observations were added to our 
sample, to obtain a better comparison between the properties of 
the different objects. 

From the collected spectra, we have compared the ratios 
of the integrated area for several relevant pairs of lines in the 
YHGs with those in AGBs and PPNe. We find that the YHGs 
studied here show, as expected, O-rich chemistry. However, the 
emission of molecules other than 13 CO in AFGL2343 is sig- 
nificantly weaker than that found for the O-rich AGBs and for 
IRC +10420. 

The comparison of the ratios of the integrated profile ar- 
eas corrected by distance has shown a clear difference between 
YHGs, whose emission is the more intense, and PPNe and AGB 
stars. Another difference, distance independent, appears when 
comparing equivalent line widths, which are much higher for 
YHGs than for the other evolved objects. At this respect, we 
note the absence of wings in the line profiles of AFGL 2343 and 
IRC +10420. 

All this supports the idea that AFGL 2343 is a YHG rather 
than a PPN, as also deduced from CO maps (Castro-Carrizo et 
al. 2007). 

We estimated the abundances of the different molecular 
species within the circumstellar envelopes of AFGL 2343 and 
IRC +10420 using the method described in Sect. 4.1. We as- 
sumed that molecular emission, apart from CO and its isotopes, 
comes from the innermost dense region found by Castro-Carrizo 
et al. (2007) and described in Table 3, which extends to several 
10 17 cm. Note that this extent is larger than that usually found 
for the same molecules in AGB stars. Dust shielding would pre- 
vent photodissociation of these molecules in the further regions 
assumed for the YHGs. On the other hand, 13 CO line emission 
in AGB CSEs and PPNe usually extends as much as 12 CO emis- 
sion; we recall that the 12 CO envelope in our sources has been 
mapped by Castro-Carrizo et al. (2007). 

The abundance found for 13 CO in AFGL 2343 is comparable 
with that usual for AGBs, while that of IRC +10420 is larger 
by a factor of 2. The abundances of other molecules found in 
IRC +10420 are comparable, in general, to those of O-rich AGB 
stars, but those of AFGL 2343 are found to be around 40 times 
lower in average. These low abundances could be due to the fact 
that these molecules are in fact abundant only in an inner region 
with a relatively low mass. 

Acknowledgements. This work has been supported by the Spanish Ministerio 
de Ciencia y Tecnologia and European FEDER funds, under grants AYA2003- 
7584 and ESP2003-04957. The contribution of AC-C was supported by the 
6fh European Community Framework Programme through a Marie-Curie Intra- 
European Fellowship. 



Humphreys, R.M., Smith, N., Davidson, K., et al. 1997, ApJ, 1 14, 2778 
Humphreys, R.M., Davidson, K. & Smith, N. 2002 ApJ, 124, 1026 
Jones, T.J., Humphreys, R.M., Gehrz, R.D., et al. 1993, ApJ, 411, 323 
Josselin, E., & Lebre, A. 2001, A&A, 367, 826 

Jourdain de Muizon, M, Cox, P., & Lequeux, J. 1990, A&AS, 83, 337 
Klochkova, V. G., Chentsov, E. L., Panchuk, V. E. MNRAS, 292, 19 
Likkel, L. 1989, ApJ 344, 350 

Lucas, R. & Guelin, M. 1990, in Submillimetre Astronomy, ed. G.D. Watt & 

A.S. Webster (Kluwer), 97. 
Maeder, A., & Meynet, G. 1988, A&AS, 76, 411 

Mauersberger, R., Guelin, M., Martin-Pintado, J., et al. 1989, A&AS, 79, 217 
Meixner, M., Ueta, J., Dayal, A., et al. 1999, ApJS, 122, 221 
Molster, F. J., Waters, L. B. F. M., Tielens, A. G. G. M.,& Barlow, M. J. 2002, 
A&A, 382, 184 

Neri, R., Kahane, C, Lucas, R., Bujarrabal, V., & Loup, C. 1998, A&AS, 130, 1 
Omont, A., Lucas, R., Morris, M., Guilloteau, S. 1993, A&A, 267, 490 
Oudmaijer, R.D., Groenewegen, M.A.T, Matthews, H.E., et al. 1996, MNRAS, 
280, 1062 

Oudmaijer, R.D. 1998, A&AS, 129, 541 

Sanchez Contreras, C, Bujarrabal, V., & Alcolea, J. 1997, A&A, 327, 689 
Skinner, C. J., Meixner, M.,& Bobrowsky, M. 1998, MNRAS, 300, 29S 
Reid, M. J., Moran, J. M., Leach, R. W., et al. 1979, ApJ, 227L, 89R 
Thevenin, F, Parthasarathy, M., Jasniewicz, G. 2000, A&A, 359, 138 



Appendix A: Spectra of other program stars. 



References 

Bachiller, R., Forveille, T., Huggins, P.J., & Cox, P. 1997a, A&A, 324, 1123 
Bachiller, R., Fuente, A., Bujarrabal, V., et al. 1997b, A&A, 319, 235 
Bujarrabal, V., Alcolea, J., Planesas, P. 1992, A&A, 257, 701B 
Bujarrabal, V., Fuente A., & Omont, A. 1994a, A&A, 285, 247 
Bujarrabal, V., Fuente A., & Omont, A. 1994b, ApJ, 421, L47 
Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., & Sanchez Contreras, C. 2001, 
A&A, 377, 868 

Castro-Carrizo, A., Lucas, R., Bujarrabal, V., Colomer, F., & Alcolea, J. 2001 
A&A, 368, L34 

Castro-Carrizo, A., Quintana-Lacaci, G, Bujarrabal, V., Neri, R., & Alcolea, J. 

2007, A&A, 465, 457 
Cernicharo, J., Alcolea, J., Baundry, A., et al. 1997, A&A, 319, 607 
de Jager, C. 1998, ARA&A 8, 145 

Hawkins, G.W., Skinner, C.J., Meixner, M., et al. 1995, ApJ, 452, 314 
Humphreys, R.M. 1991, in Wolf Rayet Stars , ed. K.A. van der Hucht & B. 
Hidayat (Kluwer), 485 



G. Quintana-Lacaci et al.: Chemistry in the CSEs around YHGs 



11 



1 

L SiO 5-4 


\ J 


i i i i i i i i 

29 SiO 2-1 i 


III 


1— I I 1 I I I I 

: 29 

- SiO 5 

f w## 


1 

-4 

1 1 1 




E 1 1 H 1 1 1 — 

7 


1 1 1 1 1 1 1 1 _ 
J CS 2-1 ': 



-500 _ 1 500 

V LS r ( km s ) 

IRC+10216 

Fig.A.2. Observed spectra towards IRC+10216. The intensity 
scale is given in units of main-beam brightness temperature (K). 




V LS r ( k m s ) 

NGC 7027 

Fig.A.3. Observed spectra towards NGC 7027. The intensity 
scale is given in units of main-beam brightness temperature (K). 



12 

6 

4 

2 



0.2 

0.1 


1.5 

1 

0.5 




G. Quintana-Lacaci et al.: Chemistry in the CSEs around YHGs 



: SiO 5-4 i 






I I I I I I I 

L 29 SiO 2-1 J 
= — I — i — i — i — — i — i — i — i — 


i i i i ii i i i 

1 \ 


- 29 

: SiO 5-4 { 

- i i i i r^Y^ 





-200 



■100 



o 1 100 



200 



V LS r ( km s ) 

RX Boo 

Fig. A.4. Observed spectra towards RX Boo. The intensity scale 
is given in units of main-beam brightness temperature (K). 

3 
2 




0.4 



i SiO 5-4 

^— i — i — i — | — i n-'^rn^ 


I - 


: 29 SiO 5-4 


i*^r j *'T J ' i — — i — i — 


—I — I — I — — I — I — I — I — 
: H 13 CN 1-0 Ii 


1 — I 1 — i 1 1 1 1 — : 

SO 2 2 — 1 1 : 



0.2 


0.04 
0.02 


-0.02 

-500 _ 1 500 

V LSR (km s ) 

TX Cam 

Fig. A.5. Observed spectra towards TX Cam. The intensity scale 
is given in units of main-beam brightness temperature (K). 



G. Quintana-Lacaci et al.: Chemistry in the CSEs around YHGs 



13 



0.05 



0.4 

0.2 


0.1 

0.05 



-0.05 
0.3 

0.2 

0.1 



0.4 

0.2 



6 

4 

2 



10 

5 





'is 

-_ CO 1-0 




■ i 1 


i i i | i . 


- 1 i 

: c 18 o 2-1 

r 1 1 1 1 1 1 


1 i I i i i 


1 ' 
L SiO 2-1 






1 i . 


B 1 1 I 

L SiO 3-2 

t , ,u 


i I, i i i i | i - 

i i 1 ' i i i i 1 i - 


i i i i i i 


SiO 5- 




i 

\ HON 1-0 


i i 

../ 


fl : 


H 1 1 1 1 M 

L HON 3-2 


H 1 1 1 1 1 Hf 

? r , , , I , - 



V LSR (km s 1 ) 



CRL 2688 



2 

1 



1.5 

1 

0.5 


1 

0.5 


2 

1.5 

1 

0.5 


2 

1.5 

1 

0.5 

3 
2 

1 


2 
1.5 

1 

0.5 


0.04 
0.02 



-0.02 
4 

2 


0.3 

0.2 

0.1 



-0.1 



H 13 CN 1-0 




SiS 15-14 



HNC 3-2 




-500 



, 
v lsr (km s ) 



500 



Fig. A.l. Observed spectra towards CRL 2688. The intensity scale is given in units of main-beam brightness temperature (K). 



